Quantum vortex dynamics in two-dimensional neutral superfluids 
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We derive an effective action for the vortex position degree-of-freedom in a superfluid by inte- 
grating out condensate phase and density fluctuation environmental modes. When the quantum 
dynamics of environmental fluctuations is neglected, we confirm the occurrence of the vortex Mag- 
nus force and obtain an expression for the vortex mass. We find that this adiabatic approximation 
is valid only when the superfiuid droplet radius R, or the typical distance between vortices, is very 
much larger than the coherence length ^. We go beyond the adiabatic approximation numerically, 
accounting for the quantum dynamics of environmental modes and capturing their dissipative cou- 
pling to condensate dynamics. For the case of an optical-lattice superfiuid we demonstrate that 
vortex motion damping can be adjusted by tuning the ratio between the tunneling energy J and 
the on-site interaction energy U. We comment on the possibility of realizing vortex-Landau-level 
physics. 
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I. INTRODUCTION 

Superfluids accommodate angular momentum by nu- 
cleating quantized vortices. ^'^ The rapid progress of re- 
search on Bose-Einstein condensed ultracold atoms in 
magnetic tiaps^iii^i^ and optical lattices-? - over the past 
decade has opened up an opportunity to study vortex 
physics in a new setting)^ i^°'^^i^^ New possibilities for 
exploring phenomena in which the quantum behavior of 
vortices is important are particularly attractive. It is 
known, for instance, that quantum fluctuations induce 
vortex lattice meltingi^ when the number of atoms per 
vortex is less than around ten, and that Boson quan- 
tum Hall statesi^ (in which vortices form incompress- 
ible quantum fluids) emerge at still lower atom densities 
(for a recent review, consult Ref. [HI)- Phase separation 
between different vortex states driven by inhomogene- 
ity in the condensate densityi^ and vortex-Peierls states 
in some optical lattices are manifestations of quantum- 
ness in vortex physics that is enhanced by geometrical 
frustration.^'' 

In spite of the long history of vortex physics in neutral 
superfluids there are still some controversies surround- 
ing the understanding of vortex quantum dynamicsJ^ 
The main concerns center on the sources of forces on 
a quantum vortex and on the effective vortex mass in 
various circumstances] "'^^'^"i^^'^'^ The consensus view(for 
completeness, however, see also Ref. based on ei- 

ther Berry phase arguments^! or hydrodynamic theory^^, 
is that a vortex has a finite mass and experiences a Mag- 
nus force proportional to the boson density in its vicinity. 
According to these theories, then, a vortex behaves like 
a massive charged particle in an effective magnetic field 
due to the background bosons. ^^'■^^ Vortex kinetic energy 
should therefore be subject to the same Landau quanti- 
zation conditions which qualitatively alter the properties 



of electrons subject to an external magnetic field. Vor- 
tex Landau quantization has nevertheless never been ob- 
served. A key aim of this work is to help point the way 
toward circumstances in which this distinctive aspect of 
vortex quantum mechanics is more likely to be exhibited. 

In early work on the properties of vortices, the vortex 
mass rriv was often treated as a purely phenomenological 
parameter More attention has been given to the micro- 
scopic underpinning of vortex mass in recent work. For 
example, Duan^^ demonstrated that the vortex mass is 
generated by the superfiuid compressibility, although the 
Magnus force was not considered. Arovas and Preire^ 
formulated a theory of vortex mass by assuming that the 
coupling between the vortex position and environment 
is dictated by gauge invariance as in electromagnetism, 
with the vortex playing the role of the charged particle, 
and phonons that of the electromagnetic field. In their 
approach the resulting minimal coupling leads to phonon 
emission by a moving vortex, which gives the vortex a 
frequency-dependent effective mass. Further progress in 
clarifying the microscopic physics of the vortex mass was 
achieved more recently in Ref. [3l] | , which approached the 
issue by using perturbation theory with a rotating pin- 
ning potential. These authors confirmed that the vortex 
mass is due to density fluctuations in the superfiuid and 
showed that the value of the mass can depend on details 
of the vortex core structure. 

In this paper, we use a different approach to derive an 
effective action for an isolated vortex. We start with 
a careful identification of the vortex-translation zero- 
modes in the Gaussian fiuctuation spectrum of a super- 
fiuid which contains a vortex. We then integrate out all 
the remaining Gaussian fluctuations, which consist of su- 
perfiuid density and phase fiuctuations distorted by the 
presence of the vortex, to obtain an effective theory which 
depends only on the vortex position degree-of-freedom. 
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In the adiabatic limit, in which the quantum dynamics of 
the environmental fluctuations can be ignored, we repro- 
duce the vortex mass predicted by Duan.— We confirm 
that the appearance of the vortex mass is due primarily 
to the coupling between vortices and condensate density 
fluctuations outside the vortex core. When quantum fluc- 
tuations of the environment modes are allowed, we find 
that the picture in which the vortex exists as a well- 
defined quantum particle with an energy-independent 
mass holds only when the vortices are in a in superfluid 
cloud with radius R , very much larger than the coherence 
length ^. In a system with many vortices the correspond- 
ing condition requires that the typical distance between 
vortices be very much larger than ^. 

We consider, in addition, the situation of a vortex in an 
optical lattice condensate and calculate the vortex mass 
and effective magnetic field numerically. Numerical re- 
sults for the vortex mass tend to agree reasonably well 
with analytical results based on a continuum model. 

Numerical results for the effective magnetic field, how- 
ever, depend strongly on the optical lattice and do not 
compare favorably with analytical results. For large vor- 
tex core sizes the discrepancy may be explained as a con- 
sequence of the large fraction of the system area that 
is occupied by vortex cores. For small vortex cores, on 
the other hand, lattice-pinning by the vortex core plays 
an important role. Taken together, these two effects im- 
ply that sharp vortex-cyclotron resonances and vortex- 



Landau-level physics are difficult to observe in optical- 
lattice Bose condensates. 

The remainder of this paper is organized as follows. In 
Sec|ll]we outline our model and present our approach to 
the evaluation of the vortex mass and effective magnetic 
field. We develop this approach analytically for the con- 
tinuum limit. In Sec. IIIII we present results of numerical 
calculations for lattice bosons. We end with conclusions 
in Sec. M\ 



II. THEORETICAL FORMULATION 

A. Effective action 

We consider a dilute interacting bosonic gas in a re- 
duced two-dimensional geometry. Our starting point is 
the partition function for the two-dimensional dilute gas 
with contact interaction (7 > 0, written as a coherent- 
state path integral over the fields 0*(x, r) and 0(x, t). It 
is given by 



Z 



-s[0*,0]/ri 



(1) 



with the Euclidean action S[(j)*^( 
grangian 



— J^^ drL and La- 



J 



2m 



.9|0(x,t)|2 



r 



(x,r) - / d: 



'(x,T)a— 0(x,r)+f[0(x) 



(2) 



The functional integration in Eq. ([T]) is over all fields pe- 
riodic on the imaginary-time axis between t — Q and 
T = fij3 = h/ksT, where is the thermal energy. 

The integral in Eq. ^ is over two-dimensional position 
X. In the context of cold atoms, the model can be real- 
ized with a sufficiently strong harmonic trapping poten- 
tial with frequency uj^ in the z-direction. The effective 
contact interaction strength g is related to the s-wave 
scattering length a of the atoms hy g = Anah'^ /{mdz), 
with m the mass of an atom. The extent of the sys- 
tem in the z-direction is denoted by dz and is related 
to the strong harmonic trapping potential in that direc- 
tion given by dz ^ ^hjmujz- We ignore the effects of 
the radial trapping potentials for now in order to avoid 
unessential complications. 

We consider field configurations close to that of a clas- 
sical vortex with counterclockwise phase winding -|-27r 
located at r„: 



Here, (/)t,(x — r„) is a local minimum of the energy func- 
tional f [(/)(x)] in Eq. ^ for any value of the vortex po- 
sition Ti,. Equivalently we can say that (^t,(x) is the sta- 
tionary solution of the Euler-Lagrange equation of the 
action in Eq. (I2|), (i.e., bS\(\)%,<l)^l 3<^* = 0) containing 
one vortex at the origin. Vortex translation in time is 
described by rt,(r). However, the fluctuating part of 
the field (5(/)(x, t) contains zero modes that also describe 
the translation of the vortex. For example, a fluctuation 
5(j} oc dx4>v describes translation of the vortex in the x- 
direction. Generally, we denote the zero modes responsi- 
ble for translation of the vortex in the x and y-direction 
by with i G {x,?/}. To avoid double-counting in 
the functional integral we have to enforce the constraint 
that the fluctuations {5(j)* , 5(j)) are orthogonal to the zero 
modes. Hence, the fluctuations are required to obey 



(^o.|'50) = / dx(/.S,(x)(50(x,r) =0 



(x, r) = 0^(x - r„(T)) -I- J0(x, r). 



(3) 



(4) 



This procedure precisely separates vortex-position and 
what we will refer to as environmental fluctuations. To 
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enforce the above constraint and, at the same time, bare vortex motion contribution Sy, an environmental- 
introduce the vortex position as a dynamical vari- fluctuation contribution S^, and a coupHng term Sc'- 
able in the function integral we use the Fadeev-Popov 
procedure^^i^. We rewrite the partition function as 

S[r, , Scf,* , (50] = Sy [r,] + Se [Scj)* , S^] + Sc [r. , Scf,* , 50] . 
Z = J d[ry]d[rW]S ((0o.|<50)) 6 {{cl^oyM) e-^['--^^*'«^l/'-> , (6) 

(5) Because £[(j)y{x — ry)] is independent of r^, the bare 

and expand the action up to quadratic order in the flue- vortex contribution to the action is due entirely to the 
tuations. The resulting actions consists of three parts, a kinetic terms: 



Sy[ry] = {(j)i,\hdr(l)u) = -[{(j}u\fldr(j)u) - {^dr(j)u\(l)u) 



ih 



Im 



9C(x) 90, (x) 



dx 



dy 



^ / drlx^iji, - y^Xi,]. (7) 



The spatial integral factor in Eq. Q, is readily evaluated for a circularly-symmetric vortex wave function 0^(x) 
/(r) exp{i6) using polar coordinates: 



/ 



dxim 



(x) 90, (x) 



dx 



dy 



drf— = 7r/^(oo) = nric 



(8) 



r 



where ric = \(f>v\'^ is the condensate density in the uniform 
superfluid far from the vortex. As a result, we obtain 



5*1, [r^] = ihimc 



dr 



drjir) 
dr 



(9) 



where e*^ is the two-dimensional Levi-Civita tensor and 
summation over repeated Cartesian indices i,j £ {x,y} 
is implied. After a Wick rotation from imaginary to 
real time (r — > it), this contribution to the action 
can be identified with the magnetic potential which ap- 
pears for charged particles in a magnetic field. When 
added to a free particle Lagrangian it implies kinetic en- 
ergy quantization and Landau level formation. The ef- 



fective magnetic length i = {hc/qBY^"^ corresponding 
to this effective magnetic potential is i^s = (27rnc)~^^^. 
The effective Aharonov-Bohm phase accumulated by a 
vortex wave function surrounding an area A is there- 
fore A/i'^ = 27tNc where Nc = Aric is the number of 
enclosed condensate particles. The correspondence be- 
tween this phase and the phase iirNy experienced by 
bosons enclosing Ny vortices is the origin of the boson- 
vortex duality properties which are usefu P^'^^ in analyz- 
ing superconductor-insulator phase transitions. 

It is convenient to express the quadratic action for 
quantum fiuctuations of the environment S'e[(S0*, <50] in 
the form 



1 f^^ f 
Se^^J dr J rfx<5*t(x) . 



OT 



n 



where the vector 5^' = [50(x, r), (50* (x, r)] and the hermitian Hamiltonian 7^ 



n = 



2m 



Ai + 2.g|0„(x)|2 



2m 



-Ai + 2g|0„(x)|^ 



(10) 



(11) 



The action contribution that describes coupling between vortex translation and and the quantum fluctuations of its 
condensate environment is 



Sc 



hp 



drdx [9i0*(x)(50(x,T) - 50*(x, T)9i0„(x)] f^(r) 
I 



(12) 



Only kinetic Berry-phase coupling terms appear here eigenfunction of the hermitian Hamiltonian Ti. 
because the vortex-translation mode is a (zero-energy) 
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In order to focus on vortex properties we next integrate 
out the environmental fluctuations Scj)* and Scj). The con- 
straints imposed by the delta-functional in the functional 
integral (Eq. ^) are conveniently applied by expanding 
the fluctuation fields in the basis of eigenfunctions of the 
Hamiltonian H: 



Xo 



(13) 



We expect 7i to have three zero-energy eigenfunctions. 
The first xop — {4'v,—<f>l)'^ is the global phase mode 
that results from the U{1) invariance of the action in 
Eq. (12). This mode makes no contribution to the dynam- 
ics of the vortex as will be shown later. The other two 
easily identifiable modes on which we focus are strictly 
speaking zero-modes only in the absence of vortex pin- 
ning, and represent vortex translation in orthogonal di- 
rections in the x ~ y plane. Hereafter, we will refer to 
these translational modes as zero modes in the discussion. 
Suitable linear combinations correspond to translations 
along translation along any direction in the x — y plane. 
We choose as basis functions Xoie{x,y} oc {di(j)y , di(j)l)'^ 
where i=x, y. One can demonstrate by direct substi- 
tution that these are indeed zero-energy eigenfunctions 
of Ti.. All other fluctuations can be expanded in terms 
of the finite-energy eigenfunctions of Ti. and corre- 
spond to density and phase fluctuations of a superfluid 
which has been distorted by the presence of a vortex. For 
the purpose of later discussions, we deflne the quantity 
Ta'a, an inner product between eigenfunctions a', a, as 
Ta'Q = ^{xli'WzlXa) , where Uz is the z-th Pauli matrix, 
with the hermitian property Ta>a = ^aa'- ^'^^ example, 
the Berry-phase coupling between the x and y compo- 
nents of the vortex displacement, considered in Eq. ([5]), 
is Txy which is purely imaginary. 



B. Adiabatic vortex dynamics and effective mass 

Now we consider vortex dynamics in the adiabatic 
approximation. This amounts to neglecting the ex- 
plicit time dependence of environmental modes in the 
action of fluctuations SeiScp* ,d(l)]. We expand the mi- 
croscopic action S in the basis of fluctuation Hamilto- 
nian Ti. eigenfunctions; that is, [S(l){x,T),S(j)* {x,t)\^ — 
J2a Ca(T)xa(x) whcrc thc Condition / dxxa{^)^ Xa' {^) = 
Sa,a' is satisfied. In the adiabatic approximation 



Se 



hp 







(14) 



The part of the action that describes the couphng be- 
comes 



Sc 



h/3 



dT 



(r)-c;(r)E,, 



Kir), (15) 



with Tj 



E i / dx [xi* (x)a.(/.. (x) - xl* (x)a,0* (x)] 

Note that the above expression is written in terms of Cq 



and c* by using the relation that S4>* (x, t) is the complex 
conjugate of <5(/)(x, t). 

Since the fluctuation action is Gaussian, we can per- 
form the integral over the environment fields c* (r) and 
Ca{T) analytically. The effective action of the vortex is 
then given by 



hp 



dT 



with the effective mass of the vortex given by 



(16) 



(17) 



Note that since the units of Via are inverse length, rn„ 
indeed has the dimension of mass. The prime on the 
summation implies omission of the zero modes, as re- 
quired by thc constraint. The vortex equation of motion 
corresponding to this action contains the Magnus force 



2Tnm, 



(18) 



where the effective magnetic field strength is related 
to the Berry-phase coupling by the expression B^s = 
2niTn[T,y\. 

The vortex mass can be estimated by approximating 
the environment modes as vortices modulated by plane- 
waves. This is a good approximation as long as we 
are interested in low energy fluctuations in which vor- 
tex motion induces slowly-varying non-uniformities. The 
phase and density environmental fluctuations modes in 
the presence of the vortex are then given respectively by 



and 



Xk 



xi 



(19) 



(20) 



where V is the area of thc system. Since these modes 
are spread over the full area of the superfluid, the cor- 
responding eigen-energies can be approximated by their 
values in the absence of a vortex: = ek = ?i^k^/2m 
and w Ek + 2gnc. When k = 0, these reduce to global 
phase and density modes respectively. 

We now use these approximations to evaluate the cou- 
pling of a vortex to its environment. Separating the con- 
densate wave function in thc presence of a vortex into 
amplitude and phase factors (0^ = e**) we have that 

dM^) = ^\<PvW' d.j{^) + e'' a.|0.(x)|, (21) 

from which it follows that thc coupling coefficients are 
given by 



ka; 



V2^ 



dx e 



-zk-x 



|<?i'«|9a;|</)«| 



(22) 
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and 



C. Beyond the adiabatic limit 



^— ^k■: 



(23) 



The vortex-core size provides the key length scale for 
vortex-condensate coupling. Since the core size corre- 
sponds to the coherence length scale ^ at which kinetic 
and interaction energy scales balance ifi^ /2m^^ = gric), 
coupling is expected to be large only when fc^ is small, 
where k = |k|. F^^ is induced by changes in wave- 
function amplitude which are non-zero only inside the 
vortex core. Choosing polar coordinates centered on the 
vortex position and approximating by ^/n^ir/^) in- 
side the vortex core and by a constant y/n^ outside the 



core with dxT — —e ^^id/dk^e ^ )/ir we obtain that 



1 



V^ricV dkx jxecorc 



dxe 



(24) 



for small fci^. For the coupling to density fluctuations 
we note that the dominant contribution comes from out- 
side the vortex cores and write dxd = — sin(0)/r as 
e^^'^{d/dkye~'^^'^)/ir'^ to obtain 



2Tmr 



X 



y/2ncV dky 7fc^ 

A/2^fc2' 



(25) 



The final form for Eq. ([25)1 is valid for fc^ small. 

Since the coupling between the vortex and its environ- 
ment must be dominated by modes with fc^ small we can 
use these results to estimate the contributions of both 
phase {6) and density [p) fluctuations to the vortex mass. 
Substituting Eq.([13]) and Eq.([l5l) into Eq. ([l7| and us- 
ing 27r/^ as an ultraviolet cut-off for the evaluation, 
we find that 



-I- « m{nc7r^'^) 



, HR/0 

32 2 



(26) 



The prefactor in this expression is the total mass of all 
boson particles inside the vortex core; the two factors 
inside square brackets are estimates of the and mP 
contributions respectively. Viewing the spatial cutoff R 
as the smaller of the system size and the distance be- 
tween vortices, we see that the vortex-mass will normally 
be dominated by the second term inside brackets, which 
represents coupling to density fluctuations, especially so 
when the fraction of the system area occupied by vortex 
cores is small. We comment later that the adiabatic ap- 
proximation provides a good description of vortex quan- 
tum mechanics only when the logarithmic enhancement 
of the vortex mass due to the its l/r phase tail plays 
a substantial role. In general both contributions to the 
vortex mass should be retained. 



We have shown that we can reproduce the commonly- 
employed vortex-particle approximation by making an 
adiabatic approximation to our effective action. The next 
step is to scrutinize the adequacy of the adiabatic approx- 
imation. We expand the fluctuating vortex and environ- 
ment fields, 



rl{T) = 



(27) 



The vortex action 



in Matsubara frequencies ijJ„ = 
is derived by integrating over Cq,,„ and c* „ as in the 
adiabatic case. The resulting action 



(28) 



is of the same form as the adiabatic approximation action 
except that the mass becomes frequency-dependent: 

■ml''{iujn) = X! -^"''"(*^»)r,.j_Q/r*fc^^ (29) 



where 



Ka^a'{it^n) = [^aSa' ,a " 2ihLL>nT a' ,a] ^ (30) 



The mass kernel captures both the renormalization of the 
effective mass at larger frequencies and the dissipative 
coupling of the vortex to the condensate environment. 
The quantity 

r„',„ = IJ dF[x^nx)xi(x) - x^'(x)xi(x)] (31) 

is the Berry-phase coupling between environment fluctu- 
ations which we have ignored in the adiabatic theory. 

Approximating, as before, density and phase fluctu- 
ations by plane waves we replace state labels a' ,a by 
momentum labels k', k. We divide the Berry-phase cou- 
pling between fluctuations Fk',k into different categories. 
Using Eqs. HH) and ^ we find that 



^ If' 



k',k 



:Sk' 



k',k, 



(32) 



and that there is no coupling among phase fluctuations 
and density fluctuations. In this approximation then, the 
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matrix Ka'^a{i'jJn) is given by Ka>^a = -K'k(«t^n)<^k',k with 
the matrix K]^ given by 

= 1,14 I k'^i) ( -^n 4" ) ■ 

To study damping, we perform a Wick rotation iujn — > 
I 



uj + iO = co^. The poles of the matrix elements iCk[a;+] 
are the Bogoliubov modes with excitation energy = 



We now derive approximate analytic results for the 
non-adiabatic vortex action. The key step is to sum over 
environment fluctuation modes to obtain an expression 
for the frequency-dependent mass kernel m{f where j G 
x,y: 



m: 



1 



2 1,2 



-^P |2 



(34) 



At small k the Bogoliubov excitation energy takes its 
phonon dispersion Hmit, Eij = hk{gnc/my^^ — Tike, 
where c is the speed of sound in the superfluid. Note 
that the final term in square brackets vanishes. 



We proceed as in the adiabatic limit, separating the 
mass into phase me and density rup fluctuation contribu- 
tions: 



me 



Eb{k) — huj — iS Eb{k) + fiLj + i5 



12 



Ebik) 



32c4 



2t4 /■2irc/C 



duj' 



2-Kc/R 



,'2 



,'2 



LU — UJ 



(35) 



Similarly, the condensate density fluctuation contribu- tion mp is 
I 



1 



1 



mir^ n 



27TC/5 



Eb{k) — hoj — i5 Eb{k) + Tilj + i5 

1 „r 1 



fi'^k^ 

^ir^kl' 



UJ — UJ 



-] + inS[uj' — uj) 



mnS^'^nc 



[2ln{2Trc/uj^) + i7T] 



(36) 



Assuming that the density-coupling is dominant we see 
that the dissipative contribution to the mass is frequency- 
independent and that the real part of the mass has a cor- 
responding logarithmic dependence on frequency. This 
result is similar to that of Arovas and Freire'^° obtained 
by different methods. The adiabatic limit result for the 



mass quoted previously corresponds to setting the fre- 
quency argument in this expression to its minimum value 
~ c/R. The Quality- factor for vortex cyclotron motion 
in the adiabatic limit is therefore 

vortex HR/O- (37) 
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Sharp cyclotron modes are expected only for very di- 
lute vortices or for isolated vortices in a large condensed 
cloud. 

We can now comment on the validity of using the adi- 
abatic approximation to describe vortex quantum fluc- 
tuations. As a criterion for the validity of the adiabatic 
theory in the continuum limit, we require that the adia- 
batic cyclotron frequency lo^ be much smaller than c/^, 
that is, 



and 



(38) 



This inequality is satisfied only if i? or the mean-distance 
between vortices is at least a thousand times ^, i.e. only if 
an extremely small fraction of the system area is occupied 
by vortex cores. This condition seems to be quite diffi- 
cult to satisfy for a cold-atom bosonic system. In the next 
section we explore the possibility of altering this compar- 
ison using the condensate tuning possibilities available in 
optical lattices. 



III. VORTEX CYCLOTRON MODES IN 
OPTICAL LATTICES 

A. Boson Hubbard Model 

Cold atoms in optical lattices offer the opportunity to 
alter vortex physics by tuning system parameters like the 
interaction strength between atoms, and the optical po- 
tential depth and period. When an optical lattice poten- 
tial is sufficiently deep, low-energy physics is described 
by a one band Hubbard model^i^ This standard Boson 
Hubbard model without an external trapping potential 
is described by the Hamiltonian 



(39) 



Here {i,j) represents summation over the nearest neigh- 
bors, J is the tunneling amplitude and U is the interac- 
tion energy between atoms on the same site, and /x is the 
chemical potential. The energy functional £ in the co- 
herent state path integral action for the Boson Hubbard 
model is correspondingly given by 



EsF — — J 



<i,j> 



+ U 



3 



3 



'3^3- 



(40) 

Mean-field configurations of the condensate are extrema 
of this energy functional and therefore satisfy the cor- 
responding Gross-Pitaevskii (GP) equation. The Boson 
Hubbard model on a 2D square lattice with lattice con- 
stant D can be viewed as a discretization of the 2D con- 
tinuum model we have been discussing heretofore with 
the correspondences 



J 



TT 9 



(42) 



Based on continuum model considerations explained ear- 
lier, we expect the adiabatic approximation to be satisfied 
when the ratio of the number of lattice sites per vortex 
to the number of lattice sites occupied by a vortex core is 
large. In searches for vortex-Landau quantization then, 
the small vortex cores produced by strong interactions 
are favored. On the other hand we know that the Bose 
condensate is destroyed by quantum fluctuations for— 
U > Uc — 4 Jzn where z is the lattice coordination num- 
ber and n is the average number of bosons per lattice 
site, (z = 4 in our case.) Using the continuum model 
correspondences we estimate that 



1)2 



U Azn 



(43) 



There should therefore be a broad range of interaction 
strengths over which superfluidity survives and vortex 
cores are concentrated in a single square lattice plaque- 
tte, even for the small values of n typical in optical lat- 
tice systems. The numerical calculations reported on be- 
low explore the possibility, suggested by the approximate 
continuum model calculations, that sharp adiabatic- limit 
vortex cyclotron motion modes might occur for U strong, 
but still smaller than Ur- 



B. Vortex Displacement Fluctuations 

In this section we describe numerical calculations based 
on the theoretical formulation outlined in Section [Hi 
which demonstrate the possibility of achieving some de- 
gree of experimental controls over the sharpness of the 
vortex cyclotron resonance signal. To find solutions of 
the square lattice Hubbard model Gross-Pitaevskii equa- 
tions which contain vortices, we impose quasi-periodic 
boundary conditions on a finite lattice with size L x L. 
We first impose a phase winding of 2tt on this area by 
using the boundary conditions 



V'(-^,y) 



-'i7r(l- 



(44) 



where the origin is chosen at the center of the com- 
putational grid which is also the center of the vortex 
core. We then minimize the energy by integrating the 
time-dependent Gross-Pitaevskii equations over imagi- 
nary time. 



n+l 



(45) 



2mL>2 



carefully restoring the normalization J d^xj-f/'f^^ P = Na 

(41) 

after each time step and starting from an ansatz in which 
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the phase winds at a constant rate around the boundary. 
Here i,j are labels for the two-dimensional lattice points, 
the label n counts steps in an imaginary time evolution 
with time-step e and Na — nL^/D^ is the number of 
atoms in the L x L computational cell. 

The energy functional for fluctuations around the vor- 
tex ground state are described by the lattice version of 
the Hermitian Hamiltonian in Eq. pT|) . Diagonalizing 
this Hamiltonian, we always find a zero-energy global 
phase mode and two other low but finite energy fiuc- 
tuation modes which we identify as vortex translational 
modes. The vortex modes do not have precisely zero- 
energy in our simulations because of lattice pinning and 
because of the boundary conditions we apply to the sim- 
ulation cell. The number of eigenvalues of Ti. is twice the 
number of atoms in the simulation cell. In order to con- 
nect approximately with an infinite system with a finite 
density of vortices, we broaden this discrete spectra into 
continua in the spectral calculations described below by 
replacing delta- functions at eigenenergies by Lorentzian 
functions with a width S ~ /2mL^. 

We report results for four different interaction 
strengths, U = 0.05J, U = 0.40J, U ^ 2.25J and 
U ~ 3.80J which we refer to below as very weak, weak, 
strong, and very strong. According to our estimates only 
the very weak and weak cases should yield vortex cores 
which extend over many square lattice plaquettes and 
this expectation is confirmed in Fig. [TJ All our calcula- 
tions were performed for a square lattice system with one 
vortex per NgUe = 46 x 46 = 2116 lattice sites. In order 
to compare with potential future optical lattice experi- 
ments we choose the number of atoms per site n = 1 in all 
the calculations described below. Our approach, which is 
designed to describe systems with weak and moderately 
strong interactions, is complementary to that taken by 
Lindner et al^ who use many-body exact diagonaliza- 
tion methods to estimate vortex properties and concen- 
trate on the case of a half-filled lattice with infinitely 
strong on-site repulsive interactions. 

The two degenerate vortex translation modes are iden- 
tified as Cx{t)S(I)ox and Cy{T)S(l)oy, corresponding to vor- 
tex displacements along the x-axis and y-axis of the lat- 
tice respectively. Here ScjjQi is a normalized eigenfunc- 
tion of the system energy functional and is obtained nu- 
merically. To make contact with the vortex coordinate 
r„ = {xv{t), yv{T)) in our continuum field theory, we can 
associate the excitation amplitudes Cx , Cy with the vortex 
coordinates , Vv using the relations 



(a) 



(b) 



CxO(PQx 
CyS(j)Oy 



-yvdy4>v 



(46) 



This identification becomes exact in the continuum limit 
where (t>y{yi — ^^{t)) = 4>v{x) — i"v • V0v for small vortex 
displacements. In this limit Cx and Cy represent precisely 
the same degrees of freedom and yy and differ only 
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FIG. 1: Condensate amplitude and phase for the vortex-state 
configuration indicated by the length and orientation of an 
arrow on a set of 12 x 12 sites near the center of the 46 x 
46 site simulation cell: a) UjJ = 0.05 b) U j J = 0.40, c) 
UjJ = 2.25 and U/J = 3.8. The continuum model estimate 
for the vortex core size in square lattice constant units is ^ = 
D{J/Uf'^ = 4.4721D, 1.5811D, 0.6667D, 0.5130D for cases 
(a), (b), (c) and (d) respectively. 



by a multiplicative constant: 



Cx 



yv 



J dx\dx<j)v 



1/2 



(47) 



In the continuum limit 



^ = / dx{f{dx9)'+(dxfr)-'^' « [7Tn,{l+\nR/0]-'/' 

Cx J 

When we turn to the Berry-phase coupling coefficients 
for the lattice model fluctuations below, they are first 
evaluated using Eq. (j3ip with the numerical eigenstates 
of the energy functional. Applying a central difference 
approximation, 9j;(/)i,(i,j) ~ [cj)^{i+l, j)-4>^{i-l, j)\/2D , 
to the vortex solutions of our lattice model are the 
lattice point indices), Eq. (|47| then allows us to appropri- 
ately renormalize those coupling coefficients that involve 
vortex-displacement coordinates. This identification of 
particular fiuctuation modes as vortex displacements is 
necessary in order to evaluate the adiabatic vortex mass 
using Eq. p7p . or its frequency-dependent non-adiabatic 
generalization using Eq. (|29p . but drops out of any di- 
rect calculation of a physical observable. In particular 
the renormalization does not influence the frequency or 
broadening of the mode which we identify as a vortex 
cyclotron resonance. Its only role is related to the physi- 
cal idcntiflcation of a particular low-frequency collective 
fluctuation coordinate as the position of a well-defined 
texture of the condensate field. The identification be- 
comes less well founded, of course, as the vortex cores 
become smaller and their locations more strongly pinned 
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FIG. 2: Comparison of the lattice- model energy functional 
eigenstate identified as a vortex displacement in the x direc- 
tion {5(j>ox (x) - circles) with the renormalized finite-difference 
derivative of the lattice-model vortex wavefunctions (solid 
line) illustrated in Fig. [1] for the same four values of U /J. 
Distance in this figure is in units of the lattice constant D of 
the optical lattice. 



by the lattice. We will return to this point frequently 
below. 

In Fig. [2l we compare the normalized vortex displace- 
ment mode (500a; (x) obtained by diagonalizing the en- 
ergy functional in Eq. [11] with the renormalized finite- 
difference approximation to the derivative, —dxff'v, of 
the vortex condensate configurations illustrated in Fig.[T] 
For the weak interaction cases (Figs. ([2^) and ([^b)) the 
close comparison between the two functions clearly jus- 
tifies the identification. In Fig. ([2^), the vortex core size 
is not very much smaller than our simulation cell size 
and we attribute differences between the two-functions to 
boundary effects. In case (Hb), the vortex core is smaller 
and boundary effects are less important. In the strong 
interaction cases, (Ht) and dUl), we see the influence of 
lattice pinning which suppresses the numerical fluctua- 
tion mode close to the vortex core, more so for stronger 
interactions. The ratio between S(j)ox{^) and —dx4>v im- 
plied by the analytical estimate, Eq. ([35]), varies from 
~ 0.31D in the weak interaction limit to ~ Q.2AD for 
the strong interaction limit and is in good qualitative 
agreement the scaling factors used to normalize —dx4>ox 
in Fig. [1 



C. Vortex Effective Mass 

As explained earlier, the vortex effective mass is due to 
coupling between the vortex displacements and their en- 
vironment which, in the Gaussian approximation we em- 
ploy here, consists of condensate phase and density fluc- 
tuations. In Fig. [3] we plot the density-of-states for con- 



FIG. 3: Density of states D(Er) of environment fluctuation 
modes as a function energy E/J for U/J — 0.05 (a), U/J = 
0.40 (b), U/J = 2.25 (c), and U/J = 3.80 (d). 



densate density and phase fluctuation modes, distorted 
by the presence of the vortex. In the absence of a vortex 
we would expect that the phase modes should have en- 
ergy — 2 J(2 — cos(fca;a) — cos(fcj,a)), varying over an 
8J range, and that the density fluctuation modes have 
energy = -I- 2U. We see in Fig. [3] that the presence 
of the vortex does not have a large influence on the over- 
all spectrum of fluctuations. The prominent features at 
E = AJ and E = 4J-f 2C/ in Fig.[3]originate from the van- 
Hove singularities of the nearest-neighbor hopping model 
on the square lattice. As the interactions strengthen, 
the phase and density fluctuation spectra separate. Note 
that the densities of states differ by a factor of two be- 
tween case (a) in which the phase and density spectra 
overlap and case (d) in which they are nearly completely 
separated. 

In Fig. [5] we plot a spectral decomposition of adiabatic- 
vortex-mass contributions dm^/dE defined by 



= / dE 



drriy 
'dE 



(49) 



In the very weak interaction case (Fig. [3I^a)) density 
and phase fluctuation contributions are not spectrally 
resolved. We see here that the mass is due mainly 
to low-energy long-wavelength fluctuation modes. For 
Figs. [D^b-d) we see separate contributions near E — 
due to phase fluctuations and near E = 2U due 
to density-fluctuations. As the interaction strengthens 
and the vortex-core radius gets smaller, the density- 
fluctuation contribution becomes relatively more impor- 
tant, in agreement with our continuum model analytic 
results. These results emphasize the fact that for moder- 
ate vortex densities, both phase and density fluctuation 
contributions need to be included in calculating the vor- 
tex mass. 
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FIG. 4: Spectral representation of adiabatic-vortex-mass con- 
tributions drriv/dE vs. mode energy E/J for U/J — 0.05 (a), 
U/J = 0.40 (b), U/J = 2.25 (c), and U/J = 3.80 (d). Where 
a separation exists, the lower-energy contribution can be asso- 
ciated with coupling to phase fluctuations and the higher en- 
ergy contribution with coupling to density fluctuations. The 
density-fluctuation contribution has a larger relative impor- 
tance when interactions are strong. 




FIG. 5: Spectral density of states of Bogoliubov excita- 
tions Db{E) as a function of E/J for interaction strengths 
U/J = 0.05 (a), U/J = 0.40 (b), U/J = 2.25 (c), and 
U/J — 3.80 (d). The Bogoliubov modes have linear disper- 
sion at long-wavelengths and a corresponding linear density 
of states at low energies. At higher energies their density 
of states shows the characteristic van-Hove singularities of 
square lattice systems with near-neighbor hopping. Linear 
dispersion holds over a wider energy range when the interac- 
tions are stronger. 



When we go beyond the adiabatic approximation, the 
kinetic couphng between environment phase and den- 
sity fluctuations is retained. This coupling leads (in 
the Gaussian fluctuation approximation) to Bogoliubov 
elementary-excitations of the environment. In the ab- 
sence of a vortex, the environment fluctuations can be 
labeled by wavevector k and the elementary excitation 
energy is given by i?_B(k) = (ek^^)^^^, with character- 
istic linear dispersion at long wavelength. For each k 
the elementary-excitation energy is therefore intermedi- 
ate between the energies of the corresponding phase and 
density modes. There is now a single branch of modes, 
instead of two-different ones. As shown in the approxi- 
mate continuum model calculations, the phase and den- 
sity contributions appear separately; their correlation 
functions have different residues at elementary-excitation 
poles. The mass becomes w-dependent and develops an 
imaginary part which captures the dissipative effect of 
environmental coupling. The sign of the contribution of 
a particular elementary excitation to the mass changes 
when CO becomes larger than its energy. 

In Fig. Owe plot the density of states Db{E) of Bogoli- 
ubov excitations in the presence of a vortex, for the usual 
set of interaction strengths. It is clear that distortion 
by the vortex has little influence on such global indica- 
tors of environment properties. For weak interactions the 
density-of-states is roughly constant except at very low 
energies, corresponding to quadratic dispersion at small 
k, whereas for stronger interactions the density-of-states 
has a linear dependence of energy, corresponding to lin- 
early dispersing phonon-like elementary excitations. 

Having characterized the environment modes, we now 
discuss the adiabatic vortex mass m„ and Berry phase 
Txy for optical lattice condensates. As explained earlier 
both of these quantities are sensitive to the conversion 
factor between the vortex coordinate Xi, and the excita- 
tion amplitude of the corresponding normalized energy 
fluctuation mode. The procedure we use to define this 
conversion factor, explained above, becomes somewhat 
arbitrary when interactions are very strong and vortices 
are strongly pinned by the lattice. The vortex cyclotron 
frequency, lUc = 2Ti\V xy\ / , is however independent of 
this conversion factor. This expression for lOc follows from 
the effective action in Eq. (I28p after Wick rotation to real 
frequencies using an adiabatic-approximation frequency- 
independent vortex mass. 

In Fig. [6l we illustrate the dependence of the adi- 
abatic vortex mass on the interaction strength U/J. 
When interactions are strong and vortex cores are small 
{U/J ^ 1), we see that the vortex mass is not much 
larger than the bare particle mass. In this strong in- 
teraction regime our numerical adiabatic vortex masses 
agree reasonably well with the continuum model expres- 
sion translated into lattice model parameters: 



my 
m 



nnJ 



7r2 IniV ln(C//J) 



32 



(50) 



where iV = 2116 is the number of lattice sites in our simu- 
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FIG. 6; Vortex mass m„ as a function of interaction strength 
U/J. The dashed-diamond line is a plot of the the analytical 
continuum model expression, Eq. (|26p . using the continuum- 
limit expression for the vortex-core size ^. The solid-circle 
line shows vortex masses evaluated numerically from Eq. (|49|l . 
The numerical calculation is dependent on the conversion fac- 
tor between the excitation amplitudes Cx{t) and Xv{t) as dis- 
cussed in the text. 



lationcell. This expression is also plotted in Fig. [51 There 
is a reasonable degree of agreement between approximate 
analytic and numerical masses in the strong interaction 
limit, even when the vortex core size becomes compara- 
ble to D, possibly because mass contributions then come 
mainly from interactions with density-fluctuations out- 
side the vortex core. At weaker interactions, where the 
vortex cores become large, the deviations from the ana- 
lytic expression are larger. In particular for very weak 
interactions the analytic expression has m,j increasing as 
, a trend that is not reproduced by our numerical cal- 
culations. We attribute this discrepancy to a violation in 
the numerical calculations of the assumption made in the 
analytic calculations that the vortex core size is small 
compared to R. At intermediate interaction strengths 
there is reasonable quantitative agreement between the 
analytic approximation, which is relatively crude when 
applied to a lattice model, and the numerical adiabatic 
mass results. 

Next we discuss numerical results for T^y^ the Berry- 
phase coupling between and Uv This coupling is the 
source of the effective magnetic field i?e// responsible 
for vortex cyclotron resonance. In the continuum model 
V^y = TTfic -Kn/D"^ ^ tt/D^. We see in Fig. [7] that the 
numerical lattice model T^y is in agreement only at inter- 
mediate values of U / J . We attribute deviations at small 
U/J to the finite size of our simulation cell, as in the 
adiabatic effective mass discussion. However, T^y also 
deviates strongly from the analytic result at large values 
of U / J . We attribute this to the fact that vortex cores in 
this regime are localized mainly in the central plaquette 
of the lattice. As a result, the condensate density vari- 
ation is almost negligible across the entire lattice. Since 



FIG. 7: Berry-phase coupling Im[ra;y] = \^xy \ between a; „ and 
function of interaction strength U /J. The dashed- 
diamond line is the data calculated from Eq. ^ for the con- 
tinuum model. The upturn at small U/J reflects the fact 
that the condensate density evaluated at the simulation cell 
boundaries is increased compared to the average condensate 
density when the vortex cores occupy a significant fraction of 
the simulation cell area. The solid-circle line is the data for 
the lattice model directly calculated from the wave function 
of the vortex displaced modes Stfiox, 54'Oy 
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FIG. 8: Vortex cyclotron frequency Ec/J as a function of in- 
teraction U /J. The dashed line shows the continuum model 
estimates while the solid line illustrates numerical lattice 
model results. 



the continuum model T^y is proportional to the difference 
between the condensate density far away from and at the 
center of the vortex core, the small values obtained for 
Txy in this limit are not surprising. The strong discrep- 
ancy compared to the continuum model in this regime is 
also due to lattice pinning effects which work against the 
identification of the low-energy fiuctuation modes with 
vortex translations. 

Before turning to the full non-adiabatic calculations, 
we first discuss the adiabatic approximation for the vor- 
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tex cyclotron frequency in the lattice model which is pro- 
portional to T^y/rriy. Using Eqs. (US]), (gl]), and (gH), the 
continuum model adiabatic cyclotron mode energy trans- 
lates into the following lattice model expression: 



8nU 



InN , ln{U/J) 
2 ~^ 2 



(51) 
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Based this expression, we should expect that the cy- 
clotron mode frequency Ed J would scale with the in- 
teraction strength nil/ J — + U/J, as indicated by the 
dashed line in Fig. [S] The discrepancies between contin- 
uum model and lattice model results illustrated in Fig. [5] 
are mainly due to discrepancy in T^y values at large U 
and due to large vortex core sizes and related discrepan- 
cies in the vortex mass at small U. At intermediate values 
of U the numerical lattice model adiabatic approximation 
and analytic continuum model adiabatic approximation 
give similar values for the vortex cyclotron frequency. 



D. Vortex Cyclotron Motion 

In Fig. [HI we show the behavior of the real and 
imaginary parts of the frequency-dependence in m^^ = 
mv{E — fiLu) for various values of [// J. The uj ^ limit 
of my{E — TiLo) is the adiabatic vortex mass discussed 
above. The contribution of particular environment Bo- 
goliubov modes to the vortex mass increases in magni- 
tude as huj approaches the energy of that mode from 
below and then changes sign at higher frequencies. The 
mass contribution vanishes for lj much larger than envi- 
ronment mode frequencies. 

The real and imaginary parts of the vortex mass are 
related by a standard Kronig-Kramers relationship. [See 
Eq. .] We see in Fig. [S] that for the weak interaction 
cases the imaginary part of the vortex-mass approaches 
a constant at small frequencies, as in the continuum case 
and that the real part of the vortex mass has a corre- 
sponding logarithmic peak at low frequencies. At strong 
interactions the long- wavelength low-energy contribution 
to Im[m^(£' — Tiuj)] is not dominant but a clear peak ap- 
pears at the energy of the van-Hove singularity in the 
Bogoliubov spectrum. This peak becomes quite strong 
for strong interactions and the real part of the vortex 
mass correspondingly becomes quite large before chang- 
ing sign near this frequency. Because the imaginary part 
of the vortex mass is relatively small and the real part is 
relatively frequency independent over a broad frequency 
regime, we conclude that the adiabatic approximation 
might be reasonable in describing the lattice-pinned vor- 
tex degree of freedom in strongly interacting optical lat- 
tices. Of course, as we have already explained, continuum 
model estimates completely fail to describe the frequency 
of the mode in this regime. 

In Fig. [ini we test for the existence of a sharp cy- 
clotron mode by examining the frequency dependence 
of the imaginary part of the vortex position-position 
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FIG. 9: Real part and imaginary part of the frequency depen- 
dent vortex mass mv{E) as a function of the energy in units 
of J for the four different interaction strengths. The irregu- 
lar frequency dependence at low frequencies in some plots is 
simply a reflection of incomplete broadening of our discrete 
Bogoliubov excitation spectrum and is without physical sig- 
nificance. 



spectral function Aq^i^] — ^^^ImG^j where G^j = 
{xi,{uj)x^{-uj)) = {yi,{uj)yi,{-uj)) . The corresponding 
adiabatic approximation in the lattice model is obtained 
by setting the off-diagonal elements between environmen- 
tal modes Ta^a' to zero in our calculations. As expected, 
the adiabatic vortex picture is reasonably accurate for the 
strongly interacting cases. For the weak interaction case 
illustrated in Fig. (9a), the vortex cyclotron mode fre- 
quency is renormalized and is strongly broadened because 
lm\my{Ey\ is large at the resonance position. The adia- 
batic cyclotron modes in cases (b) and (c) have a compli- 
cated frequency dependence because of the combined in- 
fluence of long- wavelength and van-Hove singularity envi- 
ronment modes. In Figs. (9d), the high-frequency modes 
near the van-Hove singularity are strongly dominant and 
corrections to the adiabatic approximation are weak. 



IV. CONCLUSIONS AND DISCUSSION 



In this paper we have developed an approach to vor- 
tex dynamics in neutral Bose superfluids that is based on 
careful separation of zero modes corresponding to vortex 
translation from the other fluctuations around the mean- 
field Bose condensate. Integrating out other fluctuations 
leads to a vortex action with a frequency-dependent vor- 
tex mass. In the adiabatic limit, in which frequency de- 
pendence is ignored, we recover the result of Duan^^ that 
the vortex mass diverges logarithmically with system size 
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FIG. 10: Imaginary part of vortex position-position correla- 
tion function ^oi['^] with (dashed lines) and without (solid 
hnes) the adiabatic approximation. The adiabatic approx- 
imation becomes more accurate when interactions are very 
strong (case (d)) because the lattice- weakened Berry-phase 
coupling Yxy suppresses the resonance frequency. The adia- 
batic approximation result for U/J = 3.80 coincides closely 
with its non-adiabatic counterpart. 



in a continuum model of the superfluid. The frequency 
dependence of the effective mass found in our approach is 
similar to the result obtained by Arovas and Freire^ with 
a completely different method. The frequency-dependent 
effective mass has a nonzero imaginary part which cap- 
tures the dissipation associated with vortex motion. 

In our approach Berry-phase coupling between or- 
thogonal vortex displacements, which causes the vor- 
tex to behave like a charged particle in an effective 
magnetic field, can be calculated straightforwardly. In 
the adiabatic limit the vortex undergoes undamped cy- 
clotron motion which quantum-mechanically leads to 
sharp vortex-Landau levels. Using our analytic results 
for the frequency-dependent complex effective mass we 
conclude that in the continuum limit vortex-Landau level 
physics can occur only when the system size is expo- 
nentially large compared to the vortex core size. In the 
case of many vortices, the vortex separation would also 
need to be exponentially large compared to the average 
distance between vortices. For these reasons it seems im- 
possible to achieve sharp vortex-Landau levels in an atom 
cloud unless interaction and kinetic energies are manip- 
ulated by an optical lattice potential. 

In order to explore vortex quantum mechanics in an 
optical lattice system, and as a test of our analytic con- 
tinuum model results, we have repeated our calculation 
for an optical lattice model in which we are able to inte- 
grate out condensate fluctuations numerically. Because 
of practical numerical restrictions we are able to per- 
form calculations only for system with a finite number 



(~ 2000) of lattice sites. Since the continuum model 
is expected to apply when the vortex core size is much 
larger than a lattice constant, we are not fully able to sim- 
ulate the regime of large isolated vortices for which we 
are able to obtain approximate results analytically. Nev- 
ertheless, for moderately weak interactions there is ap- 
proximate agreement between our numerical results and 
analytic estimates based on continuum model calcula- 
tions. As interactions strengthen we find that numerical 
results for the adiabatic vortex mass continue to compare 
favorably with our analytical results even when the vor- 
tex core is smaller than an optical lattice unit cell. The 
discrepancy between continuum model estimates and nu- 
merical results is much larger for the Berry-phase cou- 
pling between orthogonal vortex displacements, i.e. for 
the effective magnetic field responsible for vortex- Lan- 
dau levels. Analytic estimates of this effective magnetic 
field are in fact not easily recovered numerically. This is 
because, depending on the vortex core size, either the fi- 
nite size of the simulation cell cannot be ignored (for large 
vortex cores), or lattice pinning is important (for small 
vortices) . In the strong interaction limit we find that the 
effective magnetic field is substantially reduced, lowering 
the Landau level frequency and weakening its dissipative 
coupling to environmental modes of the condensate. Ac- 
cording to our calculations, strong interactions enhance 
the opportunity to observe well defined cyclotron mo- 
tion for this unexpected reason. Of course, the collective 
degree of freedom in this limit is not a simple vortex 
displacement because of lattice pinning effects. In ad- 
dition, because the atom density modulation associated 
with vortex motion is weaker at strong interactions, it 
may be more difficult to couple to this collective motion. 

Although our approach treats the environmental fluc- 
tuations within a Gaussian approximation, which is 
strictly speaking only valid for weakly-interacting sys- 
tems, we believe that our conclusions are more general 
on a qualitative level. This is because many of the re- 
sults we find are associated with the gapless Bogoliubov 
dispersion of the environmental modes. This gaplessness 
is a result of symmetry and is preserved for strong inter- 
actions. 

Finally, we remark that our approach is in principle 
general and can also be used to study the dissipative 
dynamics of other types of topological excitations, for 
example, skyrmions and monopoles in spinor Bose con- 
densates. 
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